Theory of the spin-torque-driven ferromagnetic resonance 
in a ferromagnet /normal- metal/ferromagnet structure 
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We present a theoretical analysis of current driven ferromagnetic resonance in a ferromagnet- 
normal-metal-ferromagnet tri-layer. This method of driving ferromagnetic resonance was recently 
realized experimentally by Tulapurkar et al. [Nature 438, 339 (2005)] and Sankey et al. [Phys. Rev. 
Lett. 96, 227601 (2006)]. The precessing magnetization rectifies the alternating current applied to 
drive the ferromagnetic resonance and leads to the generation of a dc voltage. Our analysis shows 
that a second mechanism to generate a dc voltage, rectification of spin currents emitted by the 
precessing magnetization, has a contribution to the dc voltage that is of approximately equal size 
for the thin ferromagnetic films used in the experiment. 

PACS numbers: 76.50.+g,72.25.Ba,75.75.+a,85.75.-d 



I. INTRODUCTION 

A decade ago, Slonczewskii and Berger— predicted 
that a spin-polarized current passing through a ferromag- 
net exerts a torque on its magnetic moment. The past 
decade has shown an abundance of experiments that have 
confirmed this theoretical prediction.— d£JLl Since spin- 
polarized currents are easily generated by passing an elec- 
trical current through a ferromagnet, the 'spin transfer 
torque' opens the way for all-electrical manipulation of 
nanoscale magnetic devices.— ^ 

Very recently, two groups have been able to use 
the spin torque to drive and detect ferromagnetic 
resonance in a ferromagnet-normal-metal-ferromagnet 
(FNF) trilaye n 10 ' 11 These experiments are designed such 
that the magnetization direction of one of the ferromag- 
nets is fixed by anisotropy forces, whereas the other mag- 
net is made of a softer ferromagnetic material or has 
a more symmetric shape so that its magnetization can 
more easily respond to the applied current or to an ap- 
plied magnetic field. In both experiments, an alternating 
electrical current is used to drive the ferromagnetic reso- 
nance, whereas the magnetization precession is detected 
through the dc voltage generated by rectification of the 
applied ac current by the time-dependent resistance of 
the device . 10 ' 11 The theoretical analysis of this experi- 
mental setup is the subject of this article. 

Not only does a spin-polarized current have an effect 
on the direction of the magnetization of a ferromagnet, 
a time varying magnetization also causes the flow of spin 
currents in normal metal conductors in electrical contact 
to the ferromagnet. This 'spin emission' was proposed 
by Tserkovnyak, Brataas, and Bauer as the cause of en- 
hanced damping of ferromagnetic resonance in thin fer- 
romagnetic films in good electrical contact to a normal 
metal substrate*^ It is also the mechanism underlying 
Berger's earlier prediction of the excitation of a dc volt- 
age by a precessing ma gne tization in an unbiased FNF 
trilayer— (see also Refs. Il4lll5| ). 

Spin emission affects the experiments of Refs. [1(3 and 
[Til in two different ways. First, through the enhancement 



of the damping spin emission broadens the ferromagnetic 
resonance. Second, the free layer's precessing magneti- 
zation emits alternating spin currents, which, in turn, 
generate a dc voltage through the time-varying spin- 
dependent conductance of the free layer That way, spin 
emission provides an alternative to rectification of the ap- 
plied ac current as a mechanism for the generation of a dc 
voltage in these experiments. Our calculations show that 
both consequences of spin emission appear or disappear 
together: If spin emission gives a significant contribution 
to the damping of the ferromagnetic resonance — which 
is the case for the few nm-thick free-layer ferromagnets 
used in the experiments — , then it also provides a sizable 
contribution to the measured dc voltage, and vice versa. 

In the remainder of this article we present the detailed 
theory of the electrical-current driven ferromagnetic res- 
onance needed to arrive at the above conclusion. In ad- 
dition, our theory allows us to calculate how the ferro- 
magnetic resonance frequency, the resonance width, and 
the asymmetry of the resonance lineshape are affected by 
embedding the free ferromagnetic layer into the FNF tri- 
layer. Our calculation proceeds in three parts. In Sec. ITT1 
we derive general expressions for the spin transfer torque, 
which we then apply to the calculation of the magneti- 
zation motion in Sec. IIII1 The generated dc voltage is 
calculated in Sec. IIV1 We conclude in Sec. [V] 



II. SPIN TRANSFER TORQUE 

A schematic drawing of the system we consider is 
shown in Fig. [TJ It consists of a ferromagnetic source 
reservoir, held at electric voltage V, a thin normal-metal 
spacer layer, a thin ferromagnetic layer, and a normal- 
metal drain reservoir. The direction n of the magnetiza- 
tion in the ferromagnetic source is considered to be fixed, 
whereas the direction m of the magnetization in the thin 
layer can change under the influence of an electrical cur- 
rent or an applied magnetic field. 

The nonequilibrium spin transfer torque arises from 
the discontinuity of the spin current J s across the free 
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FIG. 1: Schematic drawing of the ferromagnet-normal-metal- 
ferromagnet trilayer considered here. The left ferromagnet, 
with magnetization direction n, acts as the source reservoir. 
The right ferromagnet is the free layer. Its magnetization di- 
rection m can change in response to the applied current. The 
currents J(+) and J(— ) in the text are evaluated at the right 
and left sides of the free ferromagnetic layer, respectively. 
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-[J s (+)-J s R], 



(1) 



where </ s (+) and J s {~) are spin currents at the normal- 
metal-ferromagnet interfaces measured on the side of 
the drain reservoir and the spacer, respectively. We 
take the spacer layer and the free ferromagnet are suf- 
ficiently thin, so that all voltage drops occur across 
the ferromagnet-normal-metal interfaces, and spin relax- 
ation can be neglected^ (Note that neglecting spin re- 
laxation in the spacer layer is justified for the 10 nm 
thick Cu spacer used in the experiment of Ref. [Ill, which 
has a thickness much below the spin diffusion length in 
Cu. The spin diffusion length l s f in ferromagnets can 
be much smaller, however, and the experiments of Refs. 
Infill have a free layer thickness d comparable to Z s f , not 
d <C l s f- Still, we do not expect a strong effect of spin- 
flip scattering in this case, since the spin accumulation 
in the free layer remains fixed collinear with the direc- 
tion m of the magnetic moment, whereas the driving 
and detection of the ferromagnetic resonance depend on 
the misalignment of the two magnetic moments in the 
device*^) With these assumptions, the charge currents 
J c (±) and the spin currents J s (±) can be expressed di- 
rectly in terms of the charge and spin accumulations /z c 
and fi s in the spacer layer. For the charge and spin cur- 
rents J c (—) and J s (—) one has two sets of equations, one 
arising from the interface with the ferromagnetic source 
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Here G± = (G T ± G L )/2 and Gi + iG 2 = G n arc 
determined by the interface conductances for majority 
and minority electrons and by the mixing conductance 
for the interface between the ferromagnetic source and 
the normal-metal spacer, whereas g± = (g^ ± <7j.)/2 and 
.9i + ig% = 5T1 represent the equivalent quantities for 
the interface between spacer layer and the free ferro- 
magnet and for the interface between the free ferromag- 
net and the source. Numerical values for these conduc- 
tance coefficients have been obtained for the interfaces of 
various combinations of ferromagnetic and normal-metal 
materials^ 

The two sets of equations are slightly different because 
there are two ferromagnet-normal-metal interfaces be- 
tween the spacer layer and the drain reservoir, whereas 
there is only one interface between the spacer layer and 
the source reservoir, see Fig. [T}^ Also, in Eq. @, we 
omitted terms proportional to the time derivative n be- 
cause the magnetization of the source reservoir is held 
fixed. Similarly, for J c {+) and J s (+) we find 

J c (+) = -- [ 9+ /i c + g-fj, s ■ m] 

e 

J s (+) = ^{g-^ c + g+V s -m}m 



+ ^2 -91™ xm + 2~^ 92m - 



(4) 



Note that the charge current J c and the component J s ■ 
m of the spin current parallel to the direction of the 
magnetization of the free layer are conserved. 

The mixing conductances Gi + iGi and g\ + ig 2 de- 
scribe the coherent reflection of electrons with spin not 
collinear with the magnetization directions n and m off 
the interface with the fixed and free ferromagnetic lay- 
ers, respectively. We omitted terms that represent the 
coherent transmission of electrons with spin not collinear 
with n and m. The effect of coherent transmission is 
small for ferromagnets much thicker than the ferromag- 
netic coherence length, which is usually on the order of 
only a couple of atomic layers. We refer to Refs. [l2lfl6l for 
a theory in which these processes are included. Since the 
imaginary parts G 2 and g 2 of the mixing conductances 
are numerically small for metallic junctions (20% or less 



of Gi and g\) 
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we set G 2 and g 2 to zero in the 



and one arising from the interface with the free ferromag- 



following calculations. At the end of Sec. IIVI we discuss 
how our results are modified for finite G 2 and g 2 . 

The flow of electrical current through the FNF trilayer 
generates a spin-transfer torque only if the magnetization 
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directions n and m are not collinear. In the experiment 
of Rcfs. Uol and [III this is achieved by an applied magnetic 
field which orients the free-layer magnetization m at a 
finite angle with respect to the fixed-layer magnetization 
direction n in the absence of a current. Following Ref. 
ITT1 . we take this angle to be 90 degrees. We choose a 
right-handed set of coordinate axes (ei, 62, e^) such that 
n points along e% and m points along e-$ if no current 
is applied. The application of a current will cause m to 
deviate from e^. We'll be interested in the linear response 
regime, in which the magnetization components mi and 
m2 are proportional to the applied current J. 

With an alternating current bias, J c = J(i) = 
Re Joe Mt , Eqs. ((2]) and ((3|) give five independent equa- 
tions, from which one can solve for the five unknown 
variables, which are the charge and spin accumulations 
fi c and fi s in the spacer layer and the bias voltage V. 
Solving these to lowest order in the applied current, we 
find that two relevant components of the spin transfer 
torque |T]) are 
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(5) 
(G) 



is the free layer's ferromagnetic resonance frequency in 
the absence of electrical contact to the normal-metal 
spacer layer and the drain reservoir, whereas (lo\/oj 2 ) 
is the ratio of semi-major and semi-minor axis of the el- 
lipsoidal magnetization precession in that case. If r cq 
is dominated by the applied magnetic field H, one has 
lui = L02 = 7-ff . Rotating to the coordinate system with 
unit vectors e± and 62, the two components of r eq can 
be written 
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-ra-i (lo+ + u>- cos 4>) + toio;_ sin <f>] , 



Ten 2 — [mi(cj + — cos 0) — m2ijJ- sin (/)] , (10) 
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where lu± = (u>2 ± toi)/2 and <f>/2 is the rotation angle 
between e.\ and e\. 

With an applied ac current, J(t) = TLeJn.e luJt we can 
then solve for the magnetization components mi(t) = 
Remioe 4 "* and rri2(t) = Rem2oe luJt , with the result 



mio = mo- 



(Jo/e) (iu) + LO- sin < 



(ii) 



m 20 = m , (12) 



where we abbreviated 



G 1+ = G+ + {G 2 + -G 2 _)/ 9l . 



III. MAGNETIZATION DYNAMICS 



(7) 



The magnetization is driven out of equilibrium by the 
spin transfer torque of Eq. ((TJ. In order to solve for the 
full time-dependence of the magnetization, we use the 
Landau-Lifshitz-Gilbert equation , 25 ' 26 
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Here M is the magnetization per unit length, d is the 
thickness of the free ferromagnetic layer, 7 is the gyro- 
magnetic ratio, and a is the phenomenological Gilbert 
damping parameter. The equilibrium torque r Gq is the 
combination of the torque applied by the external mag- 
netic field and the anisotropy torque intrinsic to the ferro- 
magnet. Since we are interested in small deviations from 
equilibrium, we can expand r oq around the equilibrium 
direction m — e 3 , 
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where the frequencies uj\ and 102 are set by the energy 
cost for magnetization deviations along principal axes e[ 
and e' 2 perpendicular to e%. The constants u>\ and L02 de- 
pend on the dipolar field of the pinned layer, the demag- 
netization field and coercivity of the free layer, and the 
applied magnetic fields The geometric mean (uio^) 1 ^ 2 



f(u) = (1 + a 2 + — a^uj 2 — 2iuj(a + id + + <5_cj_ cos< 
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mo = jhG-/2dMG 1+ , 
and 
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9ll h 2 ( A _ G, 
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4de 2 M y Gi+ gi+Gi 
g!-yh 2 ( G + g x 



Ade 2 M \Gi+ gi + Gx 



(13) 
(14) 



(15) 
(16) 



The non-negative dimensionless numbers a± are the 
effective Gilbert damping parameters^ We need two 
damping parameters rather than one since the effective 
damping is anisotropic because of the presence of the sec- 
ond ferromagnet. 



IV. DC VOLTAGE 

Since we are interested in the dc voltage generated by 
the applied ac current, we need to calculate the voltage 
V(t) to second order in J(t). This implies that we need 
to solve Eqs. @ and ([3]) to first order in mi and TO2, 



V = J 



(2G X 



G A 



91 



V G i.9i+ 2giGi. 
2(gi + G 1 )g-G-m 1 
Gigigi+Gi+ 



tvm 2 G- 
2eGi+ 
eTO2«- ^ 
mo , 



(17) 
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where we abbreviated 



g 1+ = 2. 9+ + (g\ - gl)/G 1 



(18) 



The two terms in the first line of Eq. (fT"7|) . which are 
proportional to J and TO2, give an alternating contribu- 
tion to V only. The term proportional to J is the dc resis- 
tance of the device, whereas the term proportional to ra-x 
is the magnetic contribution to the admittance. (Elec- 
tronic contributions to the admittance occur at higher 
frequencies than the ferromagnetic resonance frequency 
and are not considered in our theory.) The dc voltage 
follows from the sub-leading terms in the second line of 
Eq. (|T7j) . which are proportional to Jmi and TO2TO1. The 
contribution proportional to Jm\ is rectification of the 
applied alternating current by the time-dependent con- 
ductance of the device. The contribution proportional 



to fnim\ follows from spin emission by the precessing 
magnetization of the free ferromagnet. 

The two terms contributing to the dc voltage are easily 
calculated using the results of the previous section. Using 
Eqs. pip and (fT2|) . one calculates the averages of the 
products Jmi and m^m\ over one period of the applied 
current, 



(Jmi) = 



(m 2 mi) 



m Q \Jo\ 2 
2e|/M| 2 
ml\J \ 2 uj 2 
2e 2 |/M| 2 

x [oj+ — uj- cos 



wlm/(w) + uj- sin <j> Re f(co)], 

(19) 

— (a+ + ck—)u)- sin</>]. 



The dc voltage then follows from substitution into Eq. 



v = m \J \ 2 g_G-{ gi + Gi) 
eG igi g 1+ G 1+ \f(Lo)\ 2 

x {uj 2 (2uj + a + + uj + a- + w_a_ cos0) - w_[(l + a 2 + + a + a-)uj 2 + uj 2 _ - u) 2 ^] sin^} . (20) 
In the limit a± <C 1 (which is appropriate for most experiments), Eq. (|20|) simplifies to the asymmetric Lorcntzian 

V = V f (21) 



with 



too Jo \ 2 g-G~{gi + G\) . 

V = ' 2 7=, —Fi i (2w+a+ +w+a_ + w_a_ cos 0), (22) 

4uj£eGigigi + Gi + 

I 



and 

2 22 

u = u + - 

S = a+Lu + + a_o;_ cos (f>, 

5' = 2 "°"- sin ^ . (23) 

2w+a+ + us+ct- + us-ot- cos 4> 

The asymmetry of the lineshape (|21[) depends on the 
anisotropy of the torque r oq and on the angle <j>/2 be- 
tween the principal axes and the direction n of the mag- 
netization of the fixed layer. In the experiment of Rcf. 
[Til the main contribution to r cq comes from the large 
magnetic field used to align the free layer magnetization 
perpendicular to n. This contribution is isotropic, which 
explains why no strongly asymmetric lineshapes were ob- 
served in Ref. [U The experiment of Ref . [l(] finds a sig- 
nificantly asymmetric lineshape if the applied magnetic 
field is small, the lineshapes becoming more symmetric 
at larger fields. Although this observation appears con- 
sistent with our theory, we should note that for Ref. [l(] 
the equilibrium torque T cq arising from the applied mag- 
netic field and shape anisotropy alone has <fi = and, 



hence, cannot explain an asymmetric lineshape. Ref- 
erence attributes the asymmetric lineshape to the 
imaginary part 52 of the mixing conductance which, if 
large enough, provides an alternative (but approximately 
magnetic-field independent) mechanism for an asymmet- 
ric lineshape, see the discussion below. 

The relative contributions of the rectification and the 
spin emission effects can be found by looking at the ratio 
of mo(Jmi)/e and (to2TOi)q;_, cf. Eq. (fT"7)) . For a± <C 1 
this ratio is 

toq(Jtoi) _ 2 <S - uj- sm.(j)(u> - aj )/uJa 
e(TO2TOi)a_ a ~(bj+ — uj- cos 4>) 

Since both terms in the numerator are of order 5 near the 
ferromagnetic resonance, whereas the denominator is of 
order a-OJo, the ratio (f24|) is of order 8/a-OJo. This is of 
order unity if a+ and <S_ are comparable, which happens 
precisely if the second term in Eq. (TT5)) is not small in 
comparison to the first. This, in turn, is the condition 
that spin emission gives a significant contribution to the 
total damping. Hence, we conclude that spin emission 
contributes significantly to the measured dc voltage if 
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and only if spin emission contributes significantly to the 
damping. Since (7712ml) is symmetric around u> — u)q, cf. 
Eq. (fT9| above, spin emission contributes to the symmet- 
ric part of the lineshape only. The antisymmetric part is 
due to the rectification of the applied ac current only. 

In our calculations we have neglected the imaginary 
parts g 2 and G2 of the mixing conductance because in 
metallic junctions they are known to be numerically small 
in comparison to the real parts g\ and G\. Inclusion of 
g 2 and G 2 leads to a small modification of the resonance 
frequency, because g 2 and G 2 change the gyromagnetic 
ratio 7 of the free ferromagnetic layer^ With correc- 
tions to first order in 52/51 only, the resonance frequency 
becomes 

ul = (u>% ~u?_) (25) 
4&_g 2 (ffiGi + + 2G 1 G 1+ - G.G+Y 
gi(GiG+ + g\G + - giG\+) 

More importantly, with nonzero g 2 and G%, there is a 
finite asymmetry in the lineshape even in the absence of 
magnetic anisotropy in the free layer^ 



5' 



with 



2loq [cj_ sin <j> — z + w_ cos <j))\ 
2u! + a+ + lu + &- + LU-oi- cos (p — 2za_w_ sin < 



Glg 1+ g 2 G- - 2g*G 1+ G 2 g- 
g 1 G 1 {g l + G 1 )g 1+ G- 



(26) 



(27) 



Again, our results are valid up to first order in g 2 / gi and 
G 2 /G\ only. 

We have also analyzed the case that the equilibrium 
angle between the fixed layer magnetization n and the 
free layer magnetization m is not 90 degrees. While this 
complicates the detailed expression for Vd c (w) (to the ex- 
tent that it cannot be reported here), it does not change 
our qualitative conclusions that (i) spin emission and rec- 
tification of the applied ac current have comparable con- 
tributions to the generated dc voltage if the free layer is 
thin enough that spin emission gives a sizable enhance- 
ment of the damping and (ii) the asymmetry of Vd c (w) 
around the resonance frequency lvq is small in the ratios 
w_ or g 2 1 g\ . The former ratio is small if the applied 
magnetic field is large enough to saturate the free ferro- 
magnet, whereas the latter ratio g 2 /gi is known to be 
numerically small for metallic junctions (of order 0.1 or 
less, see Refs. HHIH • 



V. CONCLUSION 

In this article we have presented a microscopic the- 
ory for the spin-torque driven ferromagnetic resonance 



in ferromagnet-normal-metal-ferromagnet trilayers. Our 
theory is inspired by the experiments of Refs. [l0 and [Til. 
In these experiments, an alternating current is used to 
drive the ferromagnetic resonance, while a generated dc 
voltage is used to detect the resonance. 

In addition to providing theoretical expressions for the 
width and asymmetry of the resonance, we are able to 
determine the relative magnitude of two physical mech- 
anisms that contribute to the dc voltage: rectification 
of the applied ac current and rectification of the spin 
currents emitted by the precessing ferromagnet. Both 
contributions are of similar magnitude for the thin fer- 
romagnetic films used in the experiments. The presence 
of two mechanisms to generate a direct response to peri- 
odic driving, rather than one, sets this class of magnetic 
devices apart from their semiconductor counterparts. 



A direct experimental probe of the two contributions 
to the dc voltage is to compare the dc voltage observed in 
spin-torque driven ferromagnetic resonance with the dc 
voltage generated in conventional magnetic-field driven 
ferromagnetic resonance in the same device. The latter 
follows from rectification of emitted spin currents only. 
Since spin emission gives a symmetric line shape around 
the resonance frequency co = u>o there should be a clear 
difference between the two methods to excite ferromag- 
netic resonance. A comparison of the magnitudes of both 
contributions would require a calibration of the ampli- 
tude at which the magnetization precesses. This can be 
achieved through a simultaneous measurement of the dc 
resistance of the device, which depends on the precession 
amplitude through the giant magnetoresistance effect. 
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